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Abstract 

A Semidefinite Programming (SDP) relaxation is an effective computational method to solve 
a Sensor Network Localization problem, which attempts to determine the locations of a group of 
sensors given the distances between some of them In this paper, we analyze and determine 
new sufficient conditions and formulations that guarantee that the SDP relaxation is exact, i.e., 
gives the correct solution. These conditions can be useful for designing sensor networks and 
managing connectivities in practice. 

Our main contribution is twofold: We present the first non-asymptotic bound on the con- 
nectivity or radio range requirement of the sensors in order to ensure the network is uniquely 
localizable. Determining this range is a key component in the design of sensor networks, and we 
provide a result that leads to a correct localization of each sensor, for any number of sensors. 
Second, we introduce a new class of graphs that can always be correctly localized by an SDP 
relaxation. Specifically, we show that adding a simple objective function to the SDP relaxation 
model will ensure that the solution is correct when applied to a triangulation graph. Since tri- 
angulation graphs are very sparse, this is informationally efficient, requiring an almost minimal 
amount of distance information. We also analyze a number objective functions for the SDP 
relaxation to solve the localization problem for a general graph. 

1 Introduction 

Graph Realization is a commonly studied topic which attempts to map the nodes in a graph G(V, E) 
to point locations in the Euclidean space based on the non-negative weights of the edges in E; that 
is, the weight of each edge corresponds to the Euclidean distance between the incident points. 
There are a number of applications of the graph realization problem [9, 13, 16,21,25]. In this paper, 
we focus on the application to Sensor Network Localization (SNL). 

A sensor network consists of a collection of sensors whose locations are unknown, and anchors 
whose locations are known. A common property of a sensor network is that each sensor detects 
others within a given connectivity (or radio) range and determines the distance from itself to these 
other sensors. Given this set of sensors and distances, the goal is to determine the exact location 
of each of the sensors. The problem becomes a graph realization problem by forming the weighted 
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undirected graph G(V,E), where the node set V represents the sensors and each non-negative 
weighted edge in E represents a known distance between two sensors. 

The SNL problem has received a lot of attention recently because of the formulation of its 
relaxation as a Semidefinite Program (SDP) [2| 111 } [24 } 130]. This formulation can find the exact 
locations of the sensors, given that the graph possesses certain properties. In this paper, present a 
number of additional sufficient conditions that guarantee unique localizability of the SDP relaxation 
of the SNL problem, i.e., conditions that ensure the SDP will give the correct solution in the space of 
a desired dimension. These conditions would be useful for designing sensor networks and managing 
connectivities in practice. 

1.1 Background 

We are given a graph G(V, E U E) in a fixed dimension d, where the nodes or points of V are 
partitioned into two sets: the set V a = {ai,...,a m } of m anchors (where m > d + 1) whose 
locations are known and the set V x = {x\, . . . , x n } of n sensors, whose locations are unknown. The 
edge set also consists of two distinct sets, E and E, where E = {(», j) : i,j G V x } is the set of edges 
between sensors, and E = {(k,j) : k G V a ,j G V x } is the set of edges between an anchor and a 
sensor. Moreover, for each G E (or (k,j) G E) the distance value between sensor i and sensor 
j (resp. anchor k and sensor j) is known as dij (resp. d^j)- The problem of finding the locations of 
the sensors can be formulated as finding points x\,X2, ■ ■ ■ ,x n G R rf that satisfy a set of quadratic 
equations: 

\\ Xi - Xj \\ 2 =d 2 ip V(i,j)eE_ 

\\ak-xjf =d 2 kj , V(k,j)eE. 1 ' 

From this, a number of fundamental questions naturally arise: Is there a localization or realization 
of Xj 's that solves this system? If there is a solution, is it unique? And is there a way to certify that 
a solution is unique? Is the network instance partially localizable, i.e., is the localization solution for 
a subset of the sensors unique? These questions were extensively studied in the graph rigidity and 
discrete geometry communities from a more combinatorial and theoretical prospective (see |15p20] 
and references therein). However, the computational aspects of these questions remained open, 
that is, the question of whether there is an efficient algorithm to numerically answer some of these 
questions was left open. 

The SDP relaxation model ([2]) and corresponding method aim to answer these questions com- 
putationally (see |11U30| ). For a G R n , the ith column of the identity matrix in R nxn , we define 
symmetric coefficient matrices := (0; — ej)(0; ej — ej) T and A^j := (a^; — ej)(cif.; —ej) T , where 
G lR, d is the vector of all zeros. The SDP relaxation can be represented as: 

maximize 

subject to Z {l . A1:d) = I d 

A ij .Z = d 2 V(i,j)€E_ (2) 

A kj .Z = 4j, V(k,j)€E 

z to. 

Here, Zn.^.j) represents the upper-left <i-dimensional principle submatrix of Z, the matrix dot- 
product refers to the sum of element-wise products A • B = AijBy, and Z >z means that the 
symmetric variable matrix Z is positive semidefinite. Note that problem ([2]) is a convex semidefinite 
program and can be approximately solved in polynomial time by interior-point algorithms. 
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One can see that the solution matrix Z G ]R,( rf +™) x of © is a matrix that can be decomposed 
into submatrices, 

The constraint Z y implies that Y y X T X. If Y = X T X, then the above formulation finds 
a matrix Z such that the columns of its submatrix X = [xi X2 ■ ■ ■ x n ] satisfy all quadratic 
equations in (P). 

Definition 1 ^4 sensor network is uniquely localizable if there is a unique X G E, dxn w/iose 
columns satisfy (pQ), and i/tere is no X £ 1R, hxn , for h > d, whose columns satisfy ([1]) and X ^ 
(X;0). In other words, there is no non-trivial extension of X G 1R, dxn into higher dimension h> d 
that also satisfies ([1]) [30]. 

Note that the notion of unique localizability is stronger than the notion of global rigidity. A 
sensor network is globally rigid only if there is a unique X G R dxn that satisfies CD), but it may 
also have a solution in a higher dimension space, that is a non-trivial extension of X G R dxn , which 
satisfies {T} [TpO]. 

The following theorem was proved in [30J : 

Theorem 1 An SNL problem instance is uniquely localizable if and only if the maximum rank 
solution of its SDP relaxation ([2]) has rank d, or equivalently, every solution matrix Z of ^ 
satisfies Y = X T X. Moreover, such a max-rank solution matrix can be computed approximately in 
polynomial time. 

The theorem asserts that the certification of a uniquely localizable network instance can be achieved 
by solving a convex optimization problem; the proof is constructive and produces a unique realiza- 
tion or localization solution for the original problem ([I]). 
The dual of the SDP relaxation ([2]) 

minimize Id*V + ^ Vijd^j + ^ Wkjdlj 
subject to (X + ^2 Vij A ij + Y w kjAkj^0 

(i,j)EE (h,j)£E 

is also useful, in that the solution to the dual tells us key properties about the primal. We define 
the dual slack matrix U G r(<*+™) *(<*+«•) as 



U = Co ()) + ^ VijAij + ^ W kjAkj- 



(ij')eS {k,j)£E 

for V G R rfxd . The dual slack matrix U is optimal if and only if it is feasible and meets the 
complementarity condition, ZU = 0. If complementarity holds, then rank(Z)+ rank(J7) < (d + n), 
and since rank(Z) > d, this means that rank({7) < n. Thus, if an optimal dual slack matrix has rank 
n, then every solution to ([2]) has rank d [30 . In fact, we have a stronger notion on localizability: 

Definition 2 A sensor network is strongly localizable if there exists an optimal dual slack matrix 
with rank n. 

Again, such a max-rank dual solution matrix can be computed approximately in polynomial time 
using SDP interior-point algorithms. 
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1.2 Our Contributions 



In this paper, we present new conditions that guarantee unique localizability of the SDP relaxation 
of the problem, i.e., conditions that ensure the SDP will give the correct solution so that the sensor 
network can be localized in polynomial time. We also enhance the relaxation such that the new 
SDP relaxation will produce a correct solution in dimension d to satisfy ([I]), even when the standard 
SDP relaxation ([2]) may not. More precisely, our result is twofold: 

1. A very popular graph in the context of sensor network localization is the unit-disk graph, 
where any two sensor points (or a sensor point and an anchor point) are connected if and 
only if their Euclidean distance is less than a given radio radius value r. It has been observed 
that when the radius (or radio range) increases, more sensors in the network can be correctly 
localized. There is an asymptotic analysis to explain this phenomenon when the sensor points 
are uniformly distributed in a unit-square [3]. In this paper, we present a non- asymptotic 
bound on the radius requirement of the points in order to ensure the network is uniquely 
localizable. Specifically, we decompose the area into sub-regions, which allows us to analyze 
the localizability of points in each sub-region, as opposed to the localizability of each indi- 
vidual point. We then determine the probability that all sensors will be localizable given a 
concentration of the sensors in a specified area. This may have practical impact by providing 
guidance on communication power ranges that ensure the network is uniquely localizable. 

2. The basic SDP localization model ([2]) is an SDP feasibility problem. An open question has 
been to determine whether adding certain objective functions to the basic model improves 
localizability of the problem; that is, the SDP feasible region contains high rank solutions, but 
with certain objectives the SDP optimal solution is guaranteed to be unique and low-rank. 
We give an affirmative answer for a generic class of graph, by identifying an objective function 
that will always result in a correct localization for this class of graphs. Our result may also 
have an influence on Compressed Sensing, which uses an objective function to produce the 
sparsest solution. Based on this idea, we present numerical results by comparing several SDP 
objective functions to illustrate their effectiveness. 

Moreover, although our theoretical analyses are based on exact distance measurements, similar 
extensions of our model (established in earlier SDP work) would be applicable to noisy distance 
data. 

1.3 Paper Organization 

The organization of this paper is as follows. First, Section [2] derives a lower bound for the con- 
nectivity radius in a sensor network that guarantees unique localizability with high probability. In 
section [3l we prove that given a triangulation (i.e., a planar, chordal and convex) graph, if the sum 
of the distances between nodes that do not have an edge between them is maximized, then the 
graph will be strongly localizable. We use this idea, and test a number of heuristic objective func- 
tions on a large number of random sensor networks to determine how well each works in practice. 
Our results for these heuristics are presented in Section [H 
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2 Bounding the Connectivity Radius 



In this section, we consider the unit-disk graph model (7JEHH] for sensor networks, where the 
Euclidean distance between any two sensor points (or a sensor point and an anchor point) is known 
(and the two points are connected) if and only if the distance between them is less than a given 
connectivity radius r. Assuming that the sensor points are randomly distributed in a region, we 
then establish a lower bound on radius r that guarantees unique localizability, with high probability, 
of the sensor network formed based on radius r. Actually, we establish a lower bound on radius r 
to ensure that the unit-disk graph is a (d+ l)-lateration graph, which would be sufficient for unique 
localizability. 

Definition 3 For some d,n > 1, the graph G(V, E) is a (d+l)-lateration graph if there exists a per- 
mutation of the points, {vr(l), 7r(2), . . . , ir(n)}, such that the edges of the sub-graph 7r(l), . . . , ir(d + 1) 
form a complete graph, and each successive point ir(j) for j > d + 2 is connected to d + 1 points 
in the set {vr(l), . . . — 1)}. This permutation of the points, it, is called a (d + l)-lateration 
ordering. 

It is shown in [M] that if a sensor network graph contains a spanning (d+l)-lateration graph and 
the points are in general position, then it is uniquely localizable. Zhu et al. [M] provide a rigorous 
proof, which is based on the intuitive concept that a complete graph of d + 1 points in general 
position can be always uniquely localized, and any point connected to d + 1 uniquely localized 
points can be also uniquely localized. 

Define r(p) to be the connectivity radius of the randomly distributed sensor points that ensures 
the network is uniquely localizable with probability at least p when the disk radius is r(p). To find 
a lower bound on r(p), we can find a connectivity radius for which the unit-disk graph G(V, E) will 
contain a spanning (d + l)-lateration graph with at least probability p. 

We approach the problem by considering a unit hypercube, 7~L = [0, l] d , which contains all the 
sensor points. We then split the region H into a grid of M equal sub-hypercubes in dimension d, 
say hi, /12, • • • , hu C H, where each sub-hypercube hi will have a volume of 1/M, and the length of 
each of its edges will be I := 1/ \[M. Without loss of generality, we can assume M = b d , where b is 
a positive integer and b > 3. Similarly, if the region considered is a hyper-rectangle in dimension d, 
we can assume M = bi ■ hi ■ ■ ■ bj, where b{ > 3 for i = 1, . . . , d are positive integers. This partition 
will allow us to analyze the localizability of sub-hypercubes in the region and their containing 
points, as opposed to analyzing the localizability of each point individually. 



2.1 Ensuring a Clique in the Graph 

Since a [d + l)-lateration ordering on the points will begin with a (d + l)-clique, we first find a 
lower bound on the radius r to ensure there exists at least one clique of d + 1 points in the graph. 

Proposition 1 Let Ti contain n points, and r > Zyfd = and M < (or equivalently 

r > -^==). Then, there exists at least one clique ofd+1 points in the unit-disk graph G(V,E). 

Proof: Note that is the length of the diagonal of each sub-hypercube hi. Thus, if r is 
lower bounded by the given value, then every point in a sub-hypercube will be connected to any 
other point in the same sub-hypercube. Furthermore, since there are at most sub-hypercubes, 
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by the pigeon-hole principle, at least one of them contains at least (d + 1) points and they must 
form a clique of d + 1 points in the unit-disk graph with given radius r. □ 

In what follows, we fixn = d ■ M + 1 = d • b d + 1. We will initialize the spanning (d + 1)- 
lateration graph construction by choosing r according to this lower bound, and let the points in the 
(d + l)-clique be the first d + 1 points in the lateration ordering. Since these points are randomly 
distributed, they must be in general position with probability one. Thus, we may assume that 
these d+1 points are anchors for the sensor network. This assumption is without loss of generality 
since our bound on the radius r established in the following sections will be much greater than the 
bound specified in Proposition [U simply because we need to ensure that not only does there exist 
a clique of d + 1 points, but also all sensor points in % form a spanning (d + l)-lateration graph 
with a high probability. 



2.2 Binomial Distribution Model 

One way to let the sensor points be randomly distributed throughout the area of T~L is to let 
the points be binomially distributed throughout each sub-hypercube of %. More specifically, the 
number of points, Yi, placed in each sub-hypercube hi, for i = 1, M, will be independently and 
binomially generated according to Y{ ~ B (n, 4?) with n = d-M+1 = d-b d +l. Once Yi is generated, 
we let these Yi sensor points be arbitrarily placed in general position within sub-hypercube hi. 

Using this binomial distribution model, let S n = Ylf=i denote the total number of points 
in the hypercube H. Since values YjS are independently and identically distributed and all sub- 
hypercubes are equally sized, the total number of points will be more or less evenly distributed in 
the entire hypercube T~L. Furthermore, by properties of the binomial distribution, 
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Thus, — > 1 almost surely and the assumption of binomially distributed sensor points through- 
out each sub-hypercube is statistically equivalent to assuming a uniform distribution of n points 
throughout the whole region T~L when M is sufficiently large. 



2.3 Connectivity Bound 

We now form further conditions on the connectivity radius r to ensure that the unit-disk graph 
G contains a spanning (d + l)-lateration graph. We have assumed that the points are binomially 
distributed in each sub-hypercube, parametrized as B (n, -p). First, r must satisfy Proposition Q3 
since it ensures a {d+ l)-clique in G. These points in the clique will represent the first d+1 points 
in the lateration ordering tt of a spanning (d + l)-lateration graph (Definition [3D. 

We construct an improved bound on the probability of localizability through an ordering of the 
hypercubes, hi G H, and hence an ordering on the points. For simplicity, we prove the following 
lemmas for the case of d = 2, and we refer to the sub-hypercubes as sub-squares. We also refer 
to (d + l)-lateration when d = 2 as trilateration. However, we note that the same analysis can be 
applied to hypercubes in higher dimensions, and our bound r > 2£\^2 in Lemmas [TH3] is analogous 
to the bound r > 2t\fd in dimension d. 
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Lemma 1 Assume that each sub-square in H £ IR 2 has at least one point, and r > U\[2. If 
the points of three sub-squares in the same row in three consecutive columns are in the trilateration 
ordering, then the points in all sub-squares in those three columns are also in the ordering. Similarly, 
if the points in three consecutive sub-squares in the same column are in the trilateration ordering, 
then the points in all sub-squares in those three columns are also in the ordering. 

Proof: First, note that the lower bound r > 2£y/2 ensures that all points in a given sub-region are 
connected to all points in a neighboring sub-region, that is, all points in a sub-region that shares 
either an edge or a point. 

For ease of explanation, consider the case that all points in the first three sub-squares in the 
first row of the grid are already in the trilateration. Since all points in the second row of the grid in 
the first three columns are within the connectivity range of all three points in the first row, these 
points are in the trilateration. Similarly, since each sub-square has at least one point, all points 
in the third row of the grid in the first three columns are within the connectivity range of three 
points in the second row, so these points are also in the trilateration. This pattern continues, until 
all points in the first three columns of the grid are in the trilateration. 

A generalization of this shows that if there are three sub-squares in the same row and in 
consecutive columns with points in the trilateration, and each sub-square has at least one point, 
then all points in the corresponding columns are also in the trilateration. 

An analogous result holds for three sub-squares in the same column and in consecutive rows. □ 

Lemma [T] states that if there are three sub-squares in a row with points in the trilateration, 
then the trilateration ordering extends to all squares in the corresponding columns. This concept 
is used below in Lemma which analyzes the cases depicted in Figure HJ 
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Figure 1: Conditions as described in Lemma [2] to ensure trilateration 



Lemma 2 Assume there is at least one point in each sub-square and r > H\[2. Then the associated 
unit- disk graph contains a spanning trilateration graph if either: 

a) There is a 3-clique in a non-corner sub-square 

b) There is a 3-clique in a corner sub-square and one of its neighbor squares has at least two points 

Proof: Again, note that r > 2lV2 ensures all points in a given sub-square are connected to 
all points in neighboring sub-squares. We show that if either of the conditions of Lemma [2] are 
satisfied, then there exists a trilateration ordering on the points in the graph. 

a) Consider the example in the left grid of Figure [H where there is a 3-clique in a non-corner 
sub-square. Let the points in this clique be the initial 3 points in the trilateration ordering. All 
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points in the sub-squares in the first three rows and first three columns are connected to this 
clique; let the points in these squares be next in the trilateration ordering. 

By Lemma [H all points in the first three rows of sub-squares are in the trilateration ordering. 
Since there are at least three columns in T~L, the same argument applies for the columns, and 
inductively, there is a trilateration ordering on the points that spreads throughout the entire 
hyperspace H. 

b) Now consider the right grid of Figure (TJ where there is a 3-clique in a corner sub-square, and 
there are at least two points in a neighboring sub-square. Let the points in this corner clique be 
the first 3 points in the trilateration ordering. All points in the three sub-squares that neighbor 
the 3-clique are connected to the clique and hence in the trilateration ordering. Next, let points 
the sub-squares in the third row and first three columns be the succeeding points in the ordering, 
followed by the sub-squares in the third column and first three rows. Using a similar argument 
as before from Lemma HJ we can construct a trilateration on the points in the graph, and all 
points are in the trilateration. 

Therefore, if the conditions of Lemma [2] hold, the associated unit-disk graph contains a spanning 
trilateration graph. □ 

The above lemma provides sufficient, but not necessary, conditions on a network for trilateration 
to exist, which implies unique localizability. Moreover, these are strict conditions for a sensor 
network, since the distribution of sensors in a network may not always ensure that there is one 
sensor in each sub-square. Thus, we extend these conditions to a more general case, and allow for the 
possibility of empty sub-squares. Clearly, too many empty sub-squares will result in a graph that 
is not uniquely localizable; also, if empty sub-squares exist, there must be restricting conditions to 
ensure the graph is not too sparse to ensure localizability. Thus, we establish additional properties 
of the graph that ensure a trilateration but allow for empty sub-squares. 

Definition 4 Two neighboring sub-squares are called adjacent neighbors if they do not share any 
edges, but share a point; neighbors that share an edge are called simple neighbors. A sub-square is 
called densely surrounded if all its simple neighbors have at least two points and one of the simple 
neighbors has at least 3 points. 



o 


o 


o 


o 


o 


o 


o 


o 


o 


o° 


o 


o 


o 




°o 


o 


o 


o 


o 


o 


o 
o 




o° 
o 


o 


o 


o° 




o° 


o 


o 


o 


o 


o 


o 

o 


o 


o 


o 


o 


s° 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 



Figure 2: Example of a grid that satisfies Lemma [3] 

Lemma 3 Assume every empty sub-square is densely surrounded and r > 2ty2. Then the asso- 
ciated unit-disk graph contains a spanning trilateration graph if there is a 3-clique in a non-corner 
sub-square. 
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Proof: Consider the grid in Figure [2J which shows an example of a non-corner 3-clique, and 
densely surrounded empty sub-squares. Notice that if an empty sub-square is densely surrounded, 
a trilateration ordering will spread around the empty sub-square: if the trilateration starts on one 
side of the empty sub-square, there are enough points on all sides of the empty sub-square for the 
ordering to continue around an empty sub-square. 

This example illustrates that empty sub-squares do not necessarily interfere with the trilatera- 
tion ordering under the given conditions. Notice that continuing the ordering among the first three 
rows of the grid, starting at the first column and adding points consecutively in the first three rows 
of the second, third and fourth columns will result in a trilateration ordering. 

Therefore, if the condition of Lemma [3] holds, the associated graph contains a spanning trilat- 
eration graph, and hence is uniquely localizable in dimension 2. □ 

We now use the fact that a sensor network containing a spanning trilateration is uniquely 
localizable [M] to establish a lower bound on the probability that the unit disk sensor network with 
radius r > U\pl is localizable. Define the two events: 

C := {There are only 3-cliques in corner sub-squares}, 
C := {There is a 3-clique in a non-corner sub-square}. 

Then, the probability that a graph with such randomly distributed points is uniquely localizable 
will be 

P {uniquely localizable} = P juniquely localizable|c| P |c| + P {uniquely localizable |C} P {C} 

> P | uniquely localizable|c}p{<5} 

Given that the total number of sub-squares is M = b 2 (for some integer b > 3), we introduce a 
parameter a := yig (or a := Jf/jS for general d) such that £ = a/y/n is the edge-length of each 
sub-square and we can use the same connectivity radius lower bound as before, now in terms of 
a, r(a) > (2a\/2) / \fn. The distribution of point number in each sub-square is binomial B On, j-A, 
and there are a total of (M — 4) non-corner sub-squares in %. Thus, the probability that there is 
a 3-clique in a non-corner sub-square is 




Let k be the number of empty sub-squares. By Lemma [21 P |uniquely localizable \k = 0, cX = 1, 
and if po = (1 — jj) n is the probability that one specific sub-square is empty, we have 

P{fc = ^}=(^)po(l-Po) M - t . 
Moreover, for any i < M — 4, we have 

/ \ M-4-i 
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From Lemma O we know 



P | uniquely localizable |/c = i, c\ > P |empty sub-squares are densely surrounded|/c = i, c| . 

The conditions of Lemma [3] require that empty sub-squares not have empty simple neighbors; thus 
we first find the probability that a sub-square does not have empty simple neighbors. Assume there 
are k empty sub-squares, say s\, S2, ■ ■ • , s&. Because of the independence assumption, these empty 
sub-squares are uniformly distributed. 

Given the empty sub-square si, the probability that S2 is not a simple neighbor of s\ is at least 
1 — j-fzr[ ) i the probability that S3 is not a simple neighbor of s\ or S2 is at least (1—2- M _ 2 )\ an d so 



on, so that the probability that no two empty sub-squares are neighbors is at least rijZi(l — jj^])- 
Moreover, the probability that an empty sub-square is densely surrounded, i.e., that all simple 
neighbors of an empty sub-square have at least two points and at least one of them has more than 
two points, is: 

p = P {All simple neighbors have at least two points} — P {All simple neighbors have exactly two points} 

= 1 - EU © (1 - irrf - [ffl (ir) 2 (1 " irrT 

Thus, the probability that all empty sub-squares are densely surrounded is 

P jempty sub-squares are densely surrounded) k = i, cj > p % ■ JJ(1 — -77^ — :)■ 

Note that the right hand side of the above equation is positive if i < M/5. Thus, we only 
consider grids with less than u := L-^/^J — 1 empty squares. Finally, we have the lower bound 
given by the following expression: 

P {uniquely localizable} 

> P I uniquely localizable \d\ P id\ 

= E^o P{uniquely localizable|/fc = i, C} P{C\k = i} P {k = i} 

— Y^i=o P{ un iquely localizable|A: = i, C} P {k = i}pq i 

— Pc o-P{^ = "} + Sr=i Pd i P {^ = *} x P |empty sub-squares are densely surrounded|fc = i, c\ 

> Pdfi P {k = 0} + yXi f ■ P 8 j P{k = i}x YTr\(i - (3) 



For different values of n (viewed as the total number of sensor points), we can find values of M, 
and thus a (where a 2 can be viewed as the average number of sensor points in each sub-square), 
such that the right hand side of Equation (|3|) is at least 0.99. Figures [3a] and I3bl show a and r 
versus the number of points n such that the right hand side of Equation [3] is at least 0.99. 

We also depicted our connectivity bound against two other bounds in Figure [H One can see 
that our and Aspnes' bounds are almost identical for any value of n. Thus, our result shows that 
the bound of Aspnes et al. in [3] (of r > 2y^g« for d = 2 ) i s true even when n is small, although 
it was initially proved to be an asymptotic bound when n is sufficiently large. Note that our bound, 
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Values of n and a that keep (3) at .99 Values of n and r that keep (3) at .99 




2.6 1 1 1 1 1 1 1 1 1 1 1 

2000 4000 6000 8000 10000 2000 4000 6000 8000 10000 

n, the number of nodes n, the number of nodes 

(a) a vs. number of nodes (b) r vs. number of nodes 



Figure 3: Bound on the Connectivity Radius 




Figure 4: Comparison of our connectivity bound to those of Javanmard |23j and Aspnes 
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while not in an analytical form, is proved for any value of n. On the other hand, we recently learned 
of another asymptotic bound that was independently developed by Javanmard and Montanari |23j. 
As seen from Figure HI their bound is actually weaker than ours and Aspnes'. 

Our connectivity result was proved for fi = [0, l] 2 , i.e., the unit square in dimension 2. The 
result can be extended to dimension d > 2. In summary, we have 

Theorem 2 Let fi £ [0, l] d be the unit hypercube in dimension d and be partitioned into a grid of 
M = b d equal sub-hypercubes, say hi, hz, ■ ■ ■ , Hm CH, where i = 1/6 is the edge length of each sub- 
hypercube. Let the number of sensor points in each sub-hypercube be independently and binomially 
generated according to B (n, jj) where n = d ■ M + 1, and let one of the sub-hypercubes contain 
d + 1 anchors. Then, if the connectivity radius satisfies r > 2iyd, the probability that the sensor 
network is uniquely localizable is given by expression 

Again, the parameter n of the binomial distribution can be viewed as the total number of sensor 
points in the region. We can also extend our result to another region T~L in dimension d into a grid 
of M equal sub-hypercubes in dimension d, say hi,h,2,.-., hu C fi, where each sub-hypercube h{ 
will have a volume of 1/M, and the length of each of its edges will be i := 1/ \[M. For example, 
we can assume M = b\ ■ 62 • • • bd, where hi > 3 for i = 1, . . . , d are positive integers. 

3 Unique Localization of Triangulation Graph 

The basic SDP localization model ([2]) is an SDP feasibility problem. When the network is not 
uniquely localizable, the max-rank of SDP feasible solutions is strictly greater than d. In practice, 
one may still be interested in finding a feasible SDP solution with rank d, representing one possible 
localization of points in R d . In this section, we show that adding an objective function that 
maximizes the sum of certain distances in a triangulation graph (in R 2 ) will produce a rank-2 SDP 
solution. The result should be applicable to d > 2. 

Definition 5 Consider a set of points V = {pi,P2 ■ ■ -Pn} E R 2 . ^4 triangulation, Tv, of the points 
in V is a subdivision of the convex hull of V into simplices (triangles) such that the edges of two 
simplices do not intersect or share a common face. 

Triangulation graphs and their properties have been studied in the literature [5^1121,126] . Bruck 
et al. [12] showed that an embedding on a unit disk graph with local angle information (angles 
between points) is NP-hard, while the same problem on a triangulation graph is not. Araujo et 
al. [5] introduced an algorithm to construct a triangulation graph from a unit disk graph with 
O (n log n) bit communications between points. 

Definition 6 For a triangulation Tv, we define a triangulation graph Gj- V (V,E) such that V = V 
and (pi,Pj) G E if and only if (pi,Pj) is an edge of a simplex in T-p. 

We formally decompose a triangulation Tv into an initial clique K3 and a set of actions 
A = {a\,a2 ■ ■ ■ a m }, where an action cii consists of adding a point and connecting it to either two 
adjacent points or two connected external points, where a point is called external if it is not strictly 
inside the convex hull of a cycle in the graph. This leads us to the following lemma, whose proof 
is omitted. 
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4 



5 



(b) The dashed blue lines are 
(a) Recursive Construction of a Triangulation the edges whose lengths will 

be maximized 



Figure 5: Triangulation Examples 



Lemma 4 A triangulation can be constructed recursively by either adding an external point that 
connects to two adjacent points of a simplex (triangle) already in Tv such that the new edges do 
not cross any existing edges (see FigureUfti, 1-4), or simply connecting two external points already 
in Tv to form a triangle (see Figure\^a, 5). 

Proof: By induction on an external point; see Figure \5&. □ 

Definition 7 In a triangulation graph, adjacent triangles are two triangles which share a common 
edge. A virtual edge exists between two points i and j when i and j belong to adjacent triangles, 
but (i,j) E. The set of virtual edges between sensors is denoted E v , and between sensors and 
anchors is denoted E v . 

Theorem 3 If we add an objective function to the SDP model ([2]) which maximizes the sum of 
the lengths of all virtual edges in a generic triangulation graph, then the rank of the optimal SDP 
solution is d and it produces the correct localization. 

Proof: We prove this via induction by assuming it is true for any triangulation graph with n 
points and then showing that it is also true for graphs with n + 1 points. It is clearly true for a 
single simplex when n = 3. 

Using the same notation as in ([2]), the primal SDP relaxation with this objective is: 



maximize 




subject to Z {1 . A1 . d) = I d 



Aij.Z = d 2 id ,y(i,j) G E 
A kj .Z = d 2 kj ,V(k,j) € E 

z y o 



(4) 
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and the dual of (|3J) is: 



minimize 



subject to V=[ q j + ^ + ^ w kj A 

^ fc J ~~ 5^ 



(5) 



(k,j)eE v 



Let X r 

number of points. By the induction assumption, the solution to ([!]) is Z r 



the optimal dual slack matrix in terms of its submatrices U r ' 



Wl dxn be the correct locations of points, where the superindex n represents the 

Id x n 
Xx n ) T (x n ) T x T 

Moreover, the optimal dual slack matrix U n satisfies U n • Z n = and has rank n; we can write 

'Tjn Tjn \ 

^ % , where U& G TR, dxd 

and U22 G R 71X ™ and >- 0. Note that the complementarity condition U n • Z n = means the 
elements of U n represent a stress on each edge such that the total force at all non-anchor points is 
zero (assuming, without loss of generality, a stress of —1 on all virtual edges). 

We decompose the triangulation graph into an initial simplex K 3 , and actions A = {01,02 • • • «m}- 
Without loss of generality, we assume the points in the first triangle are anchor points and let the 
last points added to the graph be x n+ \. For example, consider Figure E^b); let U 7 be the dual 
slack matrix on points 1-7 and assume the subgraph induced on the first 7 points is uniquely local- 
izable. When point 8 is added along with its incident edges, points (2, 4, 6, 8) form a clique (when 
including the virtual edge between 4 and 8, which is unique when its length is maximized). Thus, 
the sub-graph induced on points (2, 4, 6, 8) is uniquely localizable and its stress matrix is positive 
semi-definite with rank 1 [19J. 

Now consider the general case, where x n+ \ is the last point added to the graph. A new triangle 
is created by adding x n+ i, its adjacent triangle and the virtual edge, which forms a 4-clique. Define 
Oo to be the corresponding PSD stress matrix (with rank at least 1) on the graph formed by x n +i> 
the two points adjacent to x n +i (say, g and h) and the point with which has a virtual edge 
(say, k). We examine the case where g and h are sensors, however the case where at least one of 
them is an anchor is an easy extension. Assume f^o is normalized such that the stress associated 
with the virtual edge is —1. 



fin 



v 



Vgk "+ 
-Vgk 
-Vhk 
1 



Vhk 



Vgk 



-Vgk 
Vgh + Vg,n^ 
~Vgh 

—y g , n +i 



Vhk 



-Vhk 

-Vgh 
- Vgh + Vh,n+1 
—Vh,n+1 



1 \ 

~Vh,n+l 
-1 + Vg,n+1 + Vh,n+lJ 



Note that < (— 1 + y 9 , n +i + Vh,n+i) is easily shown and consider the updated stress matrix 
with 



TTn+l / U%2 nx i\ „ 

U22 -"VOlxn ) +U 



where £l £ K,( n+1 ) x ( n+1 ) is the stress matrix of the new edges, that is, Q(\g,h,k,n+\],\g,h,k,n+l]) = ^0 
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(a) Localization Results 



(b) Time Results 



The new matrix U n+1 will be feasible for the dual, since the elements of Qq are constructed 
specifically so that the matrix U n+1 will be of the right form, and Q y 0, U n y implies that 

u n+1 y o. 

Define 



711+1 



X n +1 

SX n ) T x n+1/ 

IT T v \ T 

y \ x n+l x n+l A -n) x n+l X n+l 



as the correct locations of the updated points. Then, 



jjn+l . Z n+1 = jjn . + ^[Oo]^(xf Xj) = 0. 

(id) 



Moreover, we can show that U22 y 0- Assume not, i.e., that there is a vector z £ R n+ such 



that 



which only holds if z T 



z T U^ 1 z = z 1 



U£ 2 




z + z T Qz = 0, 



u% 2 





z = and z T Qz = 0. Since U 2 2 y 0, this means that the first n 



elements of z are zero, i.e., 2(i :n ) = 0. Thus, 



z T Vtz = zl +1 tt n+1 = z 2 n+l (-l + y i>n+1 + yj. 



n+l) 



which implies z n+ \ = 0. Thus, z t l/J 2 +1 z = if and only if z = 0, implying U^ 1 y and 
rank(C/ n+1 ) = n + l. Therefore, rank(Z n+1 ) = d and the solution to ([4]) is exact. 
□ 



4 Heuristic Objective Function 



Based on the finding of the early section, we tested a number of heuristic objective functions for 
the DP relaxation model on a large number of random sensor networks. The following objective 
functions and strategies are used: 

1. (LSM) Initially maximize the sum of all non-edges, and do the following until either (i) the 
total error in edge lengths is sufficiently small, or {it) the maximum number of iterations has 
been reached: (a) Randomly choose an edge and minimize its length in the objective function 
(while maximizing all others) (b) If this results in a smaller total error in edge lengths, then 
minimize this edge length in the objective function; otherwise, maximize it 
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2. (IET) For each non-edge, run two SDP localizations: (a) Maximize all edge lengths, (b) 
Maximize all edge lengths, except the length of the given edge, which is minimized. If the 
second method results in a small error, then minimize the length of this non-edge in the 
objective function; otherwise, maximize it. 

3. (WEM) Run the SDP localization with no objective function to find a set of locations, X. 
Then, find the edge with maximum error. Run the SDP localizations two more times, 
with the objective functions of maximizing and minimizing the length of this edge from either 
i or j to an anchor; choose to maximize (minimize) this distance in the objective function if 
maximizing (minimizing) the distance from i or j to an anchor resulted in less error. 

4. (MAX) Maximize the sum of all the non-edge lengths 

5. (ZERO) No objective function 

The plots in Figures [6al and I6b1 show the results of each heuristic method. As the plot shows, 
(LSM) our-performed each other heuristic in terms of localization, however took much longer than 
(MAX) and (ZERO). The simple strategy of (MAX) seems to provide a good balance between the 
solution quality and time. 
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